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We present and solve a model for the vortex configuration of a disordered quantum Hall bilayer 
in the limit of strong and smooth disorder. We argue that there is a characteristic disorder strength 
below which vortices will be rare, and above which they proliferate. We predict that this can be 
observed tuning the electron density in a given sample. The ground state in the strong-disorder 
regime can be understood as an emulsion of vortex-antivortex crystals. Its signatures include a 
suppression of the spatial decay of counterflow currents. We find an increase of at least an order of 
magnitude in the length scale for this decay compared to a clean system. This provides a possible 
explanation of the apparent absence of leakage of counterflow currents through interlayer tunneling, 
even in experiments performed deep in the coherent phase where enhanced interlayer tunneling is 
observed. 



I. INTRODUCTION 

There has been much recent progress in the search 
for quantum condensed phases of quasiparticles in solids, 
such as Bose- Einstein condensates of excit ons, p olaritons, 
and magnons. A very interesting examplelEEl occurs for 
electron bilayers in the quantum Hall regime. When the 
two layers are close, and have individual filling factors 
V — 1/2, the Coulomb interactions produce a ground 
state in which electrons in one layer are correlated with 
holes in the other. The wavefunction of this state is 
that of a Bose-Einstein condensate of interlayer excitons, 
and it exhibits behaviors reminiscent of superfluidity and 
the Josephson effects: a small counterflow resistivity,!^ 
which can be understood as excitonic superfluidity, and a 
zero-bias tunneling anomaly,*^ which can be interpreted 
as a Josephson effect. However, the analogy is incom- 
plete, because neither the counterflow resistivity nor the 
width of the tunneling anomalji^l appears to vanish at 
flnite temperatures. 

Many theoretical works have suggested that these de- 
viations from conventional superfluid behaviors are con- 
nected to the presence of vortices. In a quantum Hall 
system physical and topological charges are related, so 
that random electric fields, created by the dopants, could 
induce vortices. The hypothesis that th is leads to a disor- 
dered vortex state has been used ^ l ^°m to explain features 
such as the width of the tunneling anomaly and the region 
of negative differential conductance. More recently, Fer- 
tig and collaborators have developed a strong-disorder 
model, in which t he di ssipation reflects the dynamics 
of a vortex liquid.^^^'^ Despite these potential conse- 
quences, however, there have been few attempts to pre- 
dict the vortex configuration in a bilayer. For weak, layer- 
antisymmetric disorder the appropriate model is a gauge 
glass, suggesting vortex hquids, gl asses, or con ven- 
tional superfluid states are possibilities .^^iEIElIIIEll This 
is supported by exact diagonalizatiorl^ of small systems 
with white noise disorder. 

The aim of this paper is to predict the vortex configu- 
ration of a quantum Hall bilayer, for the case of strong. 



long-range disorder, as is experimentally relevant for high 
mobility modulation doped samples. We argue that for 
a fixed disorder potential there is a characteristic value 
of the magnetic length, above which vortices prolifer- 
ate. We find that this proliferation corresponds to the 
formation of an emulsion of vortex-antivortex crystals. 
Our theory should be testable, since we estimate that 
the proliferation occurs in an experimentally accessible 
regime. Furthermore, we argue that the proliferation 
causes a dramatic suppression of the decay of counter- 
fiow currents. We find a new length scale for this de- 
cay which is one to two orders of magnitude larger than 
the corresponding length scale in the clean system. This 
provides a possible explanation of a long-standing puz- 
zle of the persistence of counterflow currents^ across an 
entire sample, in a regime where enhanced interlayer tun- 
neling conductance is observed. Such behavior is quan- 
titatively confirmed in recent experiments which show 
an area scaling for tunneling currents^^ up to the scale 
of 100/im. More generally, our work suggests that the 
quantum Hall bilayer could b e use d to study a disor- 
dered form of the "supersolid"'^^'20l that has previously 
attracted attention in superfluids, superconductors, and 
a clean bilayer model .I^S 

The remainder of this paper is structured as follows. In 
Sec.|ll]we develop a model for the vortex configuration of 
the bilayer, and identify the parameters which control the 
vortex density. In Sec. |III| we present numerical results 
for the ground state of the model, and compare these 
with a mean-field theory of an emulsion. 



IV 



In Sec. 

analyze the decay of countcrfiow currents in the ground 
state, suggest some further consequences of the emulsion, 
and discuss the role of antisymmetric disorder. Finally, 
Sec. IVl summarizes our conclusions. 



II. MODEL 

We begin by developing a model for the vortex config- 
uration, which we solve both numerically and in a mean- 
field approximation. Our starting point is the "coherence 
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network" picture/^ in which the bilayer consists of com- 
pressible puddles of electron liquid, separated by channels 
of the incompressible counterflow superfluid (see Fig. [T|). 
This is appropriate for the strong, smooth disorder pro- 
duced by dopants, which destroy s the superfluid over a 
significant fraction of the sample 1^21231 -y^g initially con- 
sider only layer-symmetric disorder, since the distance to 
the dopants is much larger than the interlayer separation. 
We focus on the simplest case of a balanced bilayer, where 
the filling fraction in each layer is = 1/2, and initially 
also neglect the small interlayer tunneling. 

This picture leads us to postulate the Hamiltonian 

ij i^i 

The first term is the electrostatic energy of an inhomoge- 
ncous charge distribution, written in terms of the charge 
Qi on the z*'* compressible puddle, and the inverse capac- 
itance matrix of the puddles = C~j^. The potential 
due to the dopants is contained in the continuous- valued 
shifts Qi which would be the optimum charges on the 
puddles in classical electrostatics. This Coulomb term 
was not considered by previous work on the coherence 
network. We will see that it is the competition between 
Coulomb energy and superfluid stiffness that controls the 
proliferation of vortices in the system. 

The second term in Eq. [T] models the energy of 
the channels. The condensate is characterized by a lo- 
cal phase 0(r), describing the interlayer phase coher- 
ence. Since this phase can wind by integer multiples of 
27r around each puddle we associate vorticities Vi with 
the puddles. The superfluid energy in the c hannels is 
H,f = J{ps/2)\Ve{r)\^d\, with stiffnes^sMI! ^ ^-i. 
As usual, Hsf leads to a vortex-vortex interaction Gij ^ 
— log Tij, and a constraint '^^Vi = 0. 

The topological defects of the condensate are merons,!^ 
which are vortices whose core corresponds to an unpaired 
electron in one layer. The meron charge is q — {e/2)(jv, 
where v is the vorticity, and cr = ±1 denotes the layer 
index of the core. Because of this relationship the two 
terms in Eq. [T] are coupled, and the charge disorder can 
drive vorticity in the channels. For those puddles with 
\q/e—[q/e]\ > 1/4 the electrostatics favors a half-electron 
charge, which is allowed if the vorticity around the puddle 
V is odd. This costs a superfluid energy proportional to 
V. Therefore, a puddle will have \v\ — 1 if this incurs 
a superfluid energy cost smaller than the electrostatic 
energy gain. 

These considerations allow us to identify the parameter 
controlling the vortex density in the percolating channels. 
For a wide range of parameters both the channel width 
Ic and puddle size dp will be on the order of the distance 
to the dopant layers, d^ « 200nm.^ (See Fig. [l]for illus- 
tration of these length scales.) The largest contribution 
to the electrostatic energy is the Coulomb interactions 
within each puddle. Thus we estimate the electrostatic 
energy gain of a vortex as E'cap ~ (l/2)(e/2)^/C, where 
C ^ dp ^ dd is the self-capacitance of the puddle. We 




FIG. 1: Schematic of a disordered quantum Hall bilayer with 
compressible puddles of electron liquid (dark-shaded areas) 
of size dp surrounded by channels of incompressible excitonic 
superfluid of size Ic- For smooth disorder dp Ic ^ dd and 
the depicted length scales are larger than the magnetic length 
k). In the limit of very strong disorder, Ic can become small 
and comparable to Iq. 



estimate the superfluid energy cost as the prefactor of 
the vortex energy, which is generally Eg ~ 27rps 1^^ . 
Thus the vortex density is controlled by the ratio 

Es dd 

Since dd is flxed by the sample, we expect the vortex 
density to vary with the magnetic length. 

In the limit of very strong disordeif^^ Ic becomes of 
the order the magnetic length Iq « 20nm, while dp re- 
mains of order dd- The vortex energy in this regime 
is Eg ~ 2nps{lc/dd), with the factor Ic/dd accounting 
for the fraction of the area occupied by the superfluid 
(up to numerical factors depending on the shapes of the 
puddles). Thus in the strong disorder limit the vortex 
density becomes independent of Iq, E^ap/Es ~ 1. We es- 
timate this numerical parameter by modeling the puddles 
as disks of radius dd ~ 200nm, and taking ps from the 
mean-field thcory^^ at zero interlayer separation. This 
gives i^cap ^ Eg^lK. 

Since our estimates of E'cap and Eg in the strong disor- 
der limit are comparable, it may be possible to vary the 
density of vortices in experiments. Decreasing Iq should 
take the system further from the (not unrealisticji-'^^'^^-') 
strong disorder limit, and so could lead to a reduction in 
the vortex density. More generally, reducing the vortex 
density requires a decrease in the capacitative energies, 
perhaps by placing gates on both sides of the sample as 
close as possible to the wells, or increasing the super- 
fiuid energy, perhaps in samples with smaller interlayer 
separation and larger tunneling. 

To predict the vortex density and configuration of the 
bilayer, we now derive and solve a Hamiltonian for the 
vorticity. For simplicity we consider Eij = 25ijEcs.p. The 
off-diagonal terms will not qualitatively affect the results. 
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because the ofF-site Coulomb interactions have a much 
shorter range than the vortex interactions Gij . The diag- 
onal elements are approximately constant, because they 
are controlled mainly by the characteristic puddle size. 
The main source of randomness is in the offset charges 

Taking e/2 as our unit of charge, we write the total 
charge on each puddle as Qi = qf^ + CT^Wi, where a = 
±1, V = 0,±1, and qf^ is the meron-free charge. In the 
ground state qf^ is the nearest even integer to qi. Thus 
the electrostatic energy of a vortex Vi on site i is 



(3) 



This is the only energy contribution which depends on 
the layer index of the core ai , and Ei can be minimized 
by setting ai = — sgTi\vi{qf' — qi)]. The distribution of 
q is broad on the scale of the charge quantization, be- 
cause the puddles contain many electrons, so that qf^ — qi 
is a uniformly distributed random variable between ±1. 
Thus, we see that the electrostatic energy takes the form 
H = Eivf, where Si varies randomly from site to site, 
with distribution P{si). In the approximation that i?cap 
is the same for all puddles, is uniformly distributed 
between ±i?cap- Note that in reality there will be some 
variation in i?cap from puddle to puddle, and the sharp 
edges in P{£i) at ±i?cap will be smoothed out. 

Combining the electrostatic and superfluid energies, we 
thus have an effective Hamiltonian for the vorticity 



G 



(4) 



We note that the random field is coupled to the pres- 
ence of vortices, independently of their sign. This differs 
from gauge glass models, where the random field couples 
directly to the vorticity. 



III. GROUND STATES 

Numerical results for the ground states of Eq. [4] are 
shown in Figs. [2]and[3] We adopt a lattice model (as in 
Ref. IT2]) where the channels are the edges of a square lat- 
tice of side L. We take the prefactor of the vortex energy 
Es to be our unit of energy. Thus Gij = V{rij) — V^(0) is 
the lattice solution to V^V{r) = — 27r(5(0), with the sin- 
gularity removed.'^ Ground states were obtained by sim- 
ulated annealing, with standard nearest-neighbor Monte 
Carlo moves. Each ground state is obtained by recording 
the lowest energy state obtained during an anneal, from 
a temperature of 0.5 to a temperature of 0.01, in steps of 
0.01. At each temperature we perform S sweeps of AL^ 
moves. As can be seen in Fig. [2| increasing 5' by a factor 
of 10'^ does not significantly change the results, so we are 
obtaining good approximations to the ground states. The 
results in Fig. |2] are quenched averages of Monte Carlo 
data obtained for different disorder realizations. 
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FIG. 2: Top panel: average ground-state vorticity \v\ as a 
function of disorder energy scale Scap, obtained by simulated 
annealing on systems of linear size L, with S Monte-Carlo 
sweeps per temperature step (see text). Each point is an av- 
erage over 25 disorder realizations. Si is taken from a uniform 
distribution of width 2i?cap. Bottom panel: corresponding 
ground state energies (crosses and hollow symbols), and in- 
teraction energy for L = 18, S — 1 (solid triangles). Solid 
curves show the corresponding results of the mean-field the- 
ory. 



From the top panel of Fig. |2] we see that the ground 
state is a uniform superfluid for small ii'cap, while vor- 
tices proliferate above a threshold E^^^. The threshold 
behavior in \v\ as a function of £^cap is sharp due to the 
discontinuity in the on-site energy distribution P{£i), and 
would in reality be rounded due to the variations of Ecap 
between puddles. 

Fig. [3] shows ground states obtained for a typical disor- 
der realization at two different strengths. These results 
show that the vortex ground states are not completely 
disordered, and are strongly suggestive of an emulsion 
of vortex-antivortex crystals. This structure appears be- 
cause the field in Eq. |4] does not dictate the sign of the 
vorticity. On the square lattice there is a minimum in 
the interaction Gq = 7r/8 = /.tc at wavevector q = (tt, tt), 
so for a uniform field < — /ic the ground state is a 
vortex-antivortex crystal.l^ Wherea s a ra ndom field cou- 
pling to Vi (as in a gauge glass modeF^n^f) competes with 
this ordering, the random field coupling to vf does not. 
It can therefore straightforwardly induce regions of the 
crystalline phase. 

The vortex density in Fig. |2] appears to be consistent 
with a mean-field theory of an emulsion. To develop such 
a theory, let us consider a mixture of two phases occu- 
pying fractions x and (1 — x) of the system, with energy 
densities fic and respectively. Without a random field 
the mean-field energy of such a mixture i^^ 



E = fif-X + Kx{l — x). 



(5) 



K is an interaction parameter, which corrects for the use 
of bulk energy densities in the first term. It will be the 
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FIG. 3: Ground states for a typical realization of the disorder 
at strengths i?cap = 0.8 (left) and 1.1 (right). Black/white 
are vorticities ±1, and gray is 0. L = 18 and S = 1000. 



only fitting parameter in the theory. 

To incorporate the random St, we interpret Eq. |5] as 
a mean-field approximation for the microscopic effective 
Hamiltonian 



J? = ^ /iXi + ^ 2Ji. 



X'iX j 



(6) 



where Xi = denotes a site in the vortex-free phase, and 
Xi = 1 one in the vortex crystal phase. The mean-field 
approximation is obtained by writing Xi ^ x + (xi — 
x), and discarding terms quadratic in the fluctuations. 
Demanding that the resulting energy agree with Eq. [5] 
allows us to relate h and the average Jij to /ic and k. 
We then incorporate the random field term from Eq. [4] 
EiXi, to obtain 



x^kN - 



2Kx)xi. (7) 



The mean-field equation is x = {xi), where () denotes an 
average in the ground state of Eq. [Tj For the uniform 
distribution of width 2i?cap for e^, we find that, when 



^^c 



K < Er 



1 - 



(8) 



and a; = otherwise. We can also compute the energy. 



N 



{Hml) = x'^K + 



E 

-2kx 2£'cap 



-dE. 



(9) 



The solid lines in Fig. |2]show the mean-field predictions 
of Eqs. [5] [8] and [9] with k = 0.4 chosen to give the 
threshold E'cap obtained numerically. As can be seen, 
this theory, with a single fitting parameter, gives a good 
account of the numerical results. Thus the ground state 
vorticity of Eq. |4] can indeed be understood in terms of 
the formation of an emulsion of vortex crystals. 



IV. DISCUSSION 

The presence of the vortex-crystal emulsion would af- 
fect counterfiow and tunneling experiments. Let us con- 



sider, in particular, the decay of a d.c. counterfiow cur- 
rent due to tunneling. Without the vortices, the super- 
fiuid phase 9 is obtained by minimizing the energy 



H = 



|V6lp - Ancos(6i) 



(10) 



where A is the tunneling strength, and n = 1/{2ttIq) is 
the electron density. A small static perturbation to the 
solution 9 = 0, such as a small counterfiow current in- 
jected at one edge, decays on the scale set by the Joseph- 
son length 



(11) 



estimatecP as 5/j,m. This means we should not ex- 
pect counterfiow currents to persist over more than a 
few microns due to leakage by interlayer tunneling (in 
other words, by the recombination of the interlayer ex- 
citons). This appears inconsistent with the experimen- 
tal observatiorpSl of an area scaling for the tunneling 
anomaly, up to length scales of 100/^m. 

With pinned vortices, we should instead consider the 
energy associated with the vorticity-free part of the 
supercurrents.'^ If we write the phase field of the vor- 
tices as 9q, we can separate out the vorticity-free phase 
field (j> = 9 — 9q. For a fixed vortex field 6*0, the ground 
state of the system is determined by a random-field XY 
model for the vorticity-free part of the system: 



|V0p - Ancos(0-l-6io) 



dh 



(12) 



This may be treated using standard techniques.'^SlSl! In 
the emulsion, the pinning phase 6*0 is disordered. We see 
that, in the limit of a vanishing correlation length for 9q, 
the tunneling field has no effect because it averages to 
zero. In our case, the vortex phase field has a correlation 
length £, ^ dd Xj, corresponding to the weak-disorder 
regime of the random field model. In this regime, the 
ground state 4> consists of domains of linear size Ldom, 
aligned with the average random field across the domain. 
The typical tunneling energy in the random field is given 
by the sum of random energies in the range ±An^^ for 
(idom/C)^ correlation areas. This gives a typical energy 
of ^{i/loY yJJjMonJCF ■ The cost in phase stiffness in 
the domain is of the order of /Os(idom)° in two dimen- 
sions. Balancing these two energies, we find the domain 
size 



L 



dom 



A. 



(13) 



We estimate that the Josephson length Aj ~ 5/im while 
the correlation length ^ ~ lOOnm. Therefore, this do- 
main size Ldom is a new length scale associated with the 
emulsion which could be one to two orders of magni- 
tudes larger than Aj in the clean system. Moreover, we 
see that static perturbations to this disordered ground 
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state ((/) ^ (/) + S(j)), such as an injected counterflow cur- 
rent, decay over this new length scale Ldom- Allowing 
for the considerable uncertainty in Aj, this decay length 
(~ O.Smm) predicted by our model is consistent with the 
apparent experimental bound^* (3> 0.1mm). This should 
be contrasted with the vortex-free state which, as men- 
tioned above, gives Aj 5/im as the decay length. 

The vortices in the emulsion will not be completely 
pinned, and hence their presence will affect the counter- 
flow superfluidity. Even if the vortices remain pinned to 
the puddles they can move a distance dp across them, 
leading to a reduction in the stiffness.!^ Thermally ac- 
tivated hopping of vortices between the puddles may 
lead to dissipation, as in previous work on the coher- 
ence network, '^■^ ^3 go that the emulsion may formally be 
a vortex liquid at finite temperatures. However, the dis- 
tribution of Bi in our model suggests a distribution of 
activation energies, in contrast to previous work. 

Direct tests of our theory may be possible in imag- 
ing experiments.!^ For example, our model predicts that 
charging lines corresponding to half-electron charges are 
common only when Eg. More generally, the iden- 

tity of physical and topological charge implies that the 
vortex configuration affects the charging spectra. 

Finally, let us revisit the role of layer-antisymmetric 
disorder. It will give additional terms in Eq. [3] which are 
proportional to ct^, leading to terms linear in Vi in the 
Coulomb gas [Eq. |4] . Provided the compressible puddles 
are effective at screening the antisymmetric disorder, the 



energy of the charge imbalance Oi will be approximately 
e!^ /Cm, where Cm dp/lo is the mutual capacitance of 
two puddles in opposite layers. This energy is a factor 
of lo/dp <C 1 smaller than i?cap, and the terms in Vi 
are small compared with those in vf. Thus while layer- 
antisymmetric disorder could affect correlations on very 
long scales, it will not affect the physical consequences 
described above, which are controlled by the scale dp. 



V. CONCLUSIONS 

In conclusion, we have developed a model of a disor- 
dered quantum Hall bilayer, in the experimentally rele- 
vant limit of strong, smooth disorder. We have argued 
that the ground state of this model can be understood 
as an emulsion of vortex-antivortex crystals. Our the- 
ory suggests that the density of the emulsion could vary 
significantly with magnetic length, and between samples, 
allowing its effects to be isolated experimentally. An im- 
portant physical consequence of the presence of such an 
emulsion (or other disordered vortex state) is a suppres- 
sion of the decay of counterflow currents, potentially ex- 
plaining the area scaling of the tunneling anomaly.^SI 
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